Silicon Nanocrystallites in SiC>2 Matrix: The Role of Disorder and Size 
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We compare, through first-principles pseudopotential calculations, the structural, electronic and 
optical properties of different size silicon nanoclusters embedded in a Si02 crystalline or amorphous 
matrix, with that of free-standing, hydrogenated and hydroxided silicon nanoclusters of correspond- 
ing size and shape. We find that the largest effect on the opto-electronic behavior is due to the 
amorphization of the embedded nanocluster. In that, the amorphization reduces the fundamen- 
tal gap while increasing the absorption strength in the visible range. Increasing the nanocluster 
size does not change substantially this picture but only leads to the reduction of the absorption 
threshold, following the quantum confinement rule. Finally, through the calculation of the optical 
absorption spectra both in a indipendent-particle and many-body approach, we show that the effect 
of local fields is crucial for describing properly the optical behavior of the crystalline case while it 
is of minor importance for amorphous systems. 

PACS numbers: 73.22-f; 71.24.+q,73.20.at; 78.67.Bf. 



I. INTRODUCTION 

It is generally accepted that the quantum confinement 
(QC), observed when electron and hole are constrained 
in one or more dimensions by a potential well, is essen- 
tial for the visible light emission in silicon nanoclusters 
(NCs) 1-3 . Nevertheless, some open questions still exist. 
For instance, it is not yet known how much the QC model 
correctly describes the dependence of the energy gap Eq 
on the NC size and how structure distortions and surface 
properties can influence the band-gap value. Again, con- 
troversial interpretations of the photoluminescence (PL) 
microscopic origin in low-dimensional Si structures still 
exist. Surface effects, in particular oxidation 4-6 , on the 
NC optical properties have been addressed in the last 
years. Many theoretical 7-9 and experimental 7 works have 
proved that interface properties have dramatic effects on 
the emission spectra of NCs. Besides, the size of the NC 
is a crucial parameter 10,11 to determine whenever such 
interface effects dominate (smaller NCs) or become neg- 
ligible (larger NCs) respect to QC. 

Various efforts have been done in order to describe optical 
absorption and emission spectra in such systems. In this 
context, inclusion of excitonic and many-body effects ap- 
pears to be of fundamental importance 12 . However, until 
now, no work has yet tackled the excitonic problem for 
this kind of systems due to the difficulty to deal with a 
large number of atoms. 

In 2003 the first theoretical ab-initio calculation con- 
sidering a small Si-NC embedded in a Si02 matrix 13 ' 14 
has been published. In this work, both the host matrix 



and the embedded nanostructure were fully relaxed. 
In the present work a more detailed first-principles calcu- 
lation of the structural, electronic and optical properties 
of Si-NC embedded in a Si02 matrix is reported. The 
two complementary cases of a perfectly crystalline and a 
completely amorphous system are taken into account. In 
particular, the effects generated by the surrounding ma- 
trix on the electronic properties of the NCs (of different 
size) are analyzed. The results are here compared with 
the corresponding ones of the free standing systems of 
same size. 

The paper is organized as follows. A description of the 
theoretical methods and of the systems is given in sec- 
tion II. The structural, electronic and optical properties 
are analyzed in section III for both crystalline (section 
III A) and amorphous (section IIIB) cases. The effects 
of the NCs size are analyzed in sec. IIIC. The excited 
state properties described through many-body methods 
are reported in section IV. Conclusions are presented in 
section V. 



II. THE METHOD 

The /3-cristobalitc (BC) SiC>2 is well known to give 
rise to one of the simplest Si/Si02 interface because of 
its diamond-like structure 15 . The crystalline structure 
has been obtained from a Sig40i28 cubic matrix (of size 
L= 14.32 A) by removing all oxygens included in a sphere 
of radius 4.4 A placed at the center of the cubic supercell. 
The results is a structure of 64 Si and 116 O atoms with 



10 Si bonded together to form a small crystalline skele- 
ton with T<j local symmetry before relaxation. In such 
core, Si atoms show a larger bond length (3.1 A) with re- 
spect to that of the Si bulk structure (2.35 A). No defects 
(dangling bonds) are present, and all the O atoms at the 
Si/Si02 interface are single bonded with the Si atoms of 
the NC. 

The optimized structure has been achieved by relax- 
ing the total volume of the cell (see Fig.l, left panel). 
This approach gives a correct description of the atomic 
average density in the NC region. Moreover, the small 
distortion induced in the host matrix (essentially due to 
the metastable nature of the BC) is practically irrele- 
vant on the NC and the interface configuration. On the 
other hand, relaxation criteria based on a fixed supercell 
volume produce, when oxygens are removed, a strong re- 
duction of the atomic average density, giving thus not 
realistic results. Together with the crystalline structure, 
the complementary case of an amorphous silica (a-Si02) 
has been considered. The glass model has been gener- 
ated using classical molecular dynamics (MD) simula- 
tions of quenching from a melt, as described in Rcf. 16 . 
The simulations have been done using semi-empirical 
ionic potentials 17 , assuming an effective quench rate of 
2.6 • 10 13 K/s. The glass match the well-connected defini- 
tion of Zachariascn 18 . The amorphous dot structure has 
been obtained starting from a Si640i28 amorphous silica 
cell by removing the 10 oxygen atoms included in a sphere 
of radius 3.5 A placed at the center of the cell, as shown 
in Fig. 1 (right panel). The relaxation of all the struc- 
tures have been performed using the SIESTA code 19 ' 20 
and Troullier-Martins pscudopotentials with non-linear 
core corrections. A cutoff of 38Ry on the plane-wave ba- 
sis set and no additional external pressure or stress were 
applied. Atomic positions and cell parameters have been 
left totally free to move. 

Figure 2 shows the radial distribution function g(r) of 
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FIG. 2: (color online) Radial distribution function (and in 
the inset its integral) for the bulk silicon (green solid line), 
the crystalline Siio (blue dashed line) and amorphous a-Siio 
(red dotted line) NCs without the matrix contribution. 



the crystalline (Siio) and amorphous (a-Siio) silicon NCs 
without considering the matrix contribution, compared 
with that of the NC made of the 8 atoms forming the 
primitive cubic cell of bulk silicon. In the inset the inte- 
grated radial distribution functions are shown. We note 
that the crystalline and the amorphous NCs present a 
nearest-neighbor distance of about 2.43 A, strained re- 
spect to the typical Si-bulk value (2.35 A), in fair agree- 
ment with the outcomes of Yilmaz et al. 21 . As expected, 
there is a good matching between the g(r) of the crys- 
talline NC and the Si-bulk case, while the long-range 
order is clearly broken in the amorphous system. 
Electronic and optical properties of the relaxed struc- 
tures have been obtained in the framework of DFT, using 
the ESPRESSO package 22 . Calculations have been per- 
formed using norm-conserving pseudopotentials within 
the LDA approximation with a Ceperley- Alder exchange- 
correlation potential, as parametrized by Perdew-Zunger. 
An energy cutoff of 60 Ry on the plane wave basis 
have been considered. To include the many-body cor- 
rections we have first calculated, starting from the DFT- 
LDA eigenvalues and eigenvectors, the GW 23 quasipar- 
ticle band structure using 10975 plane waves for the cor- 
relation part of the self-energy, and 20001 plane waves 
for the exchange part. The Bethe-Salpeter equation has 
been then solved considering an excitonic Hamiltonian 
made up of more than 65.000 transitions. The Haydock 
algorithm has been used to solve the Hamiltonian. Cal- 
culations have been performed using the EXC code 24 . 



FIG. 1: (color online) Stick and ball pictures of the final op- 
timized structure of Siio in /3-cristobalite matrix (left panel) 
or in a-Si02 glass (right panel). Red (dark gray) spheres rep- 
resent the O atoms, cyan (gray) spheres represent the Si of 
the matrix, and the yellow (gray) thick sticks represent the Si 
atoms of the NC. 



III. INDIPENDENT- PARTICLE RESULTS 

In parallel to the crystalline and amorphous Siio/Si02 
systems, three other structures have been studied: (i) 
the pure Si02 matrix (in the same phases as used 
for the embedded NC calculations), (ii) the isolated 



2 



NCs as extracted from the relaxed NC-silica complexes 
and capped by hydrogen atoms (Siio-H), (hi) the NCs 
together with the first interface oxygens extracted as 
in point (ii), and then passivated by hydrogen atoms 
(Siio-OH). In the last two cases only the hydrogen atoms 
have been relaxed. The goal is to distinguish between 
the properties that depend only on the embedded NC 
from those that are instead influenced by the presence 
of the matrix. The comparison of the results relative to 
different passivation regimes (H or OH groups) could 
give some insight on the role played by the interface 
region. 



A. The crystalline cluster in /3-cristobalite matrix 

From the analysis of the relaxed Siio/Si02 supercell 
emerges that the NC has a strained structure with respect 
to bulk Si 21 ' 25 ' 26 (see Fig. 2), while the /3-cristobalite ma- 
trix is strongly distorted near the NC and reduces pro- 
gressively its stress far away from the interface 27 . The 
electronic and optical properties of the relaxed structures 
of Siio/Si02, Si02 bulk, Siio-H and Siio-OH have been 
computed. Figure 3 shows the comparison between the 
band structures (energy levels) of all these systems. We 
see, first of all, a strong reduction of the Siio/Si02 energy 
gap (E G =1.77 eV) with respect to both bulk Si0 2 (5.44 
cV) and the isolated NC passivated by H atoms (4.66 eV); 
finally, as confirmed by Refs. 4 > 5 > 7 , the passivation by OH 
groups tends to red-shift the energy spectrum {Eq= 1.6 




FIG. 3: Band structure along high symmetry points of the 
BZ for the crystalline Siio NC in /3-cristobalite (a) , compared 
with the band structure of /3-cristobalite bulk (b), the energy 
levels for the isolated Siio NC passivated by H (c) and by OH 
groups (d). The arrows show the HOMO-LUMO band gap 
values. The units are in eV. 



cV). The smaller gap with respect to the Si02 bulk case 
is clearly due to the formation of confined states within 
the NC, lying between the Si02 band edges, as displayed 
in Fig. 3a and 3b. These states are fundamentally dif- 
ferent from those of the free-standing hydrogenated NC 
(Fig. 3c). In the embedded NC states, indeed, we find 
a strong contribution due to the presence of interfacial 
oxygen atoms. This is also confirmed by the resemblance 




FIG. 4: (color online) Kohn-Sham orbitals at 10% of their 
max. amplitude for the Siio crystalline embedded NC. Top 
panel: HOMO state is represented in blue (gray), HOMO-l 
in azure (light gray), HOMO-2 in lightest blue (lightest gray). 
Bottom panel: LUMO is represented in dark blue (dark gray) . 
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FIG. 5: (color online) DFT-RPA imaginary part of the di- 
electric function for the Siio crystalline NC in /3-cristobalitc 
matrix (solid line) compared with that of isolated Siio-H NC 
(long-dashed line), the beta-cristobalite bulk (short-dashed 
line), and the Siio-OH NC (dotted line). The inset shows an 
enlargement of the spectra at low energies. 





Si0 2 


Siio/SiOa 


Siio-OH 


Siio-H 


Crystalline 


5.44 


1.77 


1.60 


4.66 


Amorphous 


5.40 


1.41 


1.55 


1.87 



TABLE I: DFT HOMO-LUMO gap values (in eV) for the 
crystalline and amorphous silica, embedded, OH-terminated, 
H-terminated Siio NC. 



the transitions can not be ascribed neither to the matrix 
region nor to the pure NC. In this region the peaks are 
originated by the interplay between the embedded NC 
and the matrix; in particular they are due to the presence 
of the oxygen atoms located at the interface region. This 
is confirmed by the similarity that exists, in this energy 
region, between the Siio/Si02 and the Siio-OH spectra, 
as depicted in the insertion of Fig. 5. The presence of 
OH groups thus induces a complex structure in the part 
of the spectrum below 4 eV that can be compared with 
the peaks of the Siin/Si02 structures. 



existing between the gap values of the embedded NC and 
of the NC simply passivated by the OH groups (Fig. 3d). 
Eg is thus almost completely determined by the barrier 
provided by the first shell of oxygen atoms. The typical 
behavior of the bulk dispersed states is still recognizable 
far away from the band edges, where the NC influence 
tends to vanish. The square modulus contour plots of 
the highest occupied molecular orbital (HOMO) and the 
two states just below it (HOMO-1, HOMO-2), for the 
Siio crystalline NC in the BC matrix are reported in Fig. 
4. The supercell relaxation breaks the symmetry of the 
system, and consequently the degeneracy of the three 
states at the top of the valence band. Therefore these 
states present slightly different eigenvalues but preserve 
the character of the original symmetry. In accordance 
with Ref. 25 , these states are mainly localized on part 
of the NC, in particular on the Si-Si bonds and on the 
interface oxygens, while the lowest unoccupied molecu- 
lar orbital (LUMO) is localized on the interface and over 
all the entire NC region. It is clearly shown in Fig. 4 
that the first shell of oxygens around the NC are able to 
completely trap the NC density charge, thus forming a 
strong barrier that is responsible for the QC effect. 

In Figure 5 the imaginary part of the dielectric func- 
tions (calculated within the DFT-RPA approach, using 
the dipolc approximation and neglecting the local-field 
(LF) effects) of Sii /SiO 2 , Sii -H and Sii -OH NCs, and 
of the BC bulk are reported. A detailed analysis of the re- 
sults allows to identify three distinct zones in the spectra. 
The first one, above 6 eV, where the Siin/Si02 and BC 
bulk are very similar. The second one, between ~4 eV 
and ~6 eV, where a clear contribution due to the embed- 
ded Si-NC is present. The third one, below 4 eV, where 



B. The amorphous cluster in a glass 

The NC, when formed in the glass, completely loses 
memory of the starting tetrahedral symmetry (as shown 
in the previously discussed Fig. 2). 

No dangling bonds are present at the NC surface while 
some bridge-bonded oxygens appear (that are not present 
in the crystalline case). Despite the fact that the dra- 




FIG. 6: Band structure along high symmetry points of the 
BZ for the a-Siio/a-Si02 system (a) compared with the band 
structure of the glass bulk (b), the energy levels for the iso- 
lated a-Siio NC passivated by H (c) and by OH groups (d). 
The arrows show the HOMO-LUMO band gap values. The 
units are in eV. 
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matic structural changes with respect to the crystalline 
case, similar considerations can be done concerning the 
optoelectronic properties. While the amorphization pro- 
cess does not affect the behavior of the glass matrix, it 
strongly reduces the energy gap of the isolated a-Siio-H 
NC, of the a-Siio-OH NC, and of the composite a-Siio/a- 
Si0 2 system (see Fig. 6). The HOMO-LUMO gaps are 
summarized in table I. 

We note that, contrary to the crystalline case, the amor- 
phous hydroxided NC shows a gap value slightly greater 




FIG. 7: (color online) Kohn-Sham orbitals at 10% of their 
max. amplitude for the a-Siio amorphous embedded NC. Top 
panel: HOMO state is represented in red (gray), HOMO-1 in 
orange (light gray), and HOMO-2 in white. Bottom panel: 
LUMO is represented in dark red (dark gray). 



than when embedded. Also in this case the states at 
the valence and conduction band edges are localized at 
the interface; in particular they stem from the OH pas- 
sivation. Obviously the reduction of the local symmetry 
induces a splitting in the energy levels and, as a conse- 
quence, a more uniform distribution of states. 
The HOMO, HOMO-1, HOMO-2 and LUMO states (see 
Fig. 7) follow the shape of the NC, resulting strongly de- 
formed with respect to the crystalline case, but are still 
clearly confined by the interface oxygen shell. 
In the calculated optical absorption spectra (see Fig. 8) 
the three regions due to the hosting matrix, Si-NC, and 
interface, are still clearly distinguishable. 

The most important differences between the calculated 
optical absorption spectra for both the amorphous and 
the crystalline embedded NCs are found in the energy 
range between 2 and 4 eV. In this region the amorphous 
system shows more intense peaks with respect to the crys- 
talline case, suggesting a possible higher emission in the 
visible range. 

Despite the fact that the possibility to form amorphous 
NC in silica has been recently explored 28 , experimen- 
tal measurements on single Si-NCs with diameters of the 
order of 1 nm have not yet been performed. Thus, a 
straightforward comparison of our results with experi- 
mental data is not possible. Besides, the comparison 
with other works 29 sustains the idea that the strong de- 
formation of the NC is determinant for the absorption 
strength at low energies. This idea is also supported by 
the fact that, for larger systems, when the shape of the 
NCs tends to be spherical and the distortion is usually 
less pronounced 11,30 , crystalline and amorphous systems 
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FIG. 8: (color online) DFT-RPA imaginary part of the dielec- 
tric function for the a-Siio NC in a glass matrix (solid line) 
compared with that of isolated, hydrogenated a-Siio-H NC 
(long-dashed line), the glass matrix (short-dashed line), and 
the hydroxided a-Siio-OH NC (dotted line). The inset shows 
an enlargement of the spectra at low energies. 
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produce more similar absorption spectra. This topic will 
be discussed in the next subsection. 



C. Size Effects 

From a theoretical point of view, the simplest model 
for the QC is provided by the particle-in-a-box scheme, 
in which the box size is given by the NC diameter and 
the potential barrier represents the host insulating ma- 
trix. When QC effect dominates over other quantum 
phenomena, for an infinite barrier potential, we have 
Ea{R) = ec+A/R Q , where e G is the bulk-silicon band 
gap, R is the radius of the NC, A is a positive constant 
and a < 2 31 . This model, however, is often too simple 
to describe correctly the experimental results. 
Experimentally, several factors contribute to make the 
interpretation of measurements of the optical gaps a dif- 
ficult task. For instance, samples show a certain distribu- 
tion in the NC size, that is difficult to be determined. In 
this case it is possible that the observed PL peak does not 
correspond exactly to the mean size but instead to the 
largest PL rate. Again, NC synthesized by using different 
techniques often show different properties in size, shape 
and in the interface structure. Finally, also the NC-NC 
interactions can play a significant role on the emission 
spectrum. 

Recent measurements on samples with diameters in the 
range of 2-10 nm have reported 32-38 a PL peak that is 
blue-shifted by AE G (R) ~ R~ a , with 1 < a < 1.5. 
Thus, the model described above is far from being able 
to reproduce the measurements. In order to improve it, 
the finitcness of the potential barrier and also other cor- 
rections should be taken into account. For example, the 
inclusion of the electron-hole Coulomb interaction, that 
contributes a further dependence in blue-shifting the 
energy, results in a more gradual dependence of the op- 




FIG. 9: (color online) Stick and ball pictures of the final opti- 
mized structure of Si32 NC in a /3-cristobalite matrix (left) and 
a-Si32 NC in a a-Si02 glass (right). Red (dark gray) spheres 
represent the O atoms, cyan (gray) spheres represent the Si 
of the matrix, and the yellow (gray) thick sticks represent the 
Si atoms of the NC. 



tical gap on the dot radius, as observed experimentally. 
As mentioned above, our investigation explores the small 
dot size limit, with a diameter range within 1 nm. In 
this range, most of the NC atoms are positioned at the 
interface, where the effects of stress and oxidation are 
stronger. For bigger NCs, when Si-bulk states emerge 
and the surfacc-to-volume ratio decreases, these effects 
are limited; in this case we expect a response that is inde- 
pendent of the specific geometrical configuration. Thus, 
in order to understand all the interface-related phenom- 
ena, an investigation of the small size limit becomes nec- 
essary. 

To model NC of increasing size, we enlarge the hosting 
matrix so that the separation between NCs is still around 
1 nm, that we consider enough to correctly describe the 
stress localized around each NC and to avoid an overlap- 
ping of the NC states that can be caused by the applica- 
tion of the periodic boundary conditions 13 . These con- 
siderations lead us to choose a BC matrix of 648 atoms 
(in a cubic cell with lattice parameter of 2 nm) with a 
core of 32 silicon atoms, having an average diameter of 
1 nm. The relaxed Si32 NC shows a strained interface 
with average Si-Si bond-lengths of 2.6 A, while the sil- 
icon atoms in the core have a less strained bond-length 
of 2.43 A. Note that Sii /SiO 2 and Si 32 /Si0 2 have ap- 
proximatively the same ratio between the NC and the 
matrix volumes. The amorphous embedded a-Si3 2 NC 
is again produced by MD annealing and successive ab- 
initio relaxation. The final structures are shown in Fig. 
9. Once more, bridge-bonded oxygens are present at the 
interface and we find that the number of Si-O-Si bonds 
increases with the dimension of the NC (8 bonds in this 
case with respect to the 3 present in the smaller NC) in 
nice agreement with previous results obtained by differ- 
ent methods 11,39 . 

E G values calculated in DFT for Si 32 /Si0 2 , Si 32 -OH 
and Si3 2 -H structures are reported in Tab. II. We can 
note that, while Eq of Si3 2 -H is smaller with respect to 
Siio-H due to the QC, for the other structures it strongly 
depends on the interface configuration. As a conse- 
quence, for both the Si3 2 /Si0 2 and Si3 2 -OH systems, 
Eq increases when the phase is crystalline, and decreases 
when the corresponding structures are amorphous. We 
note that, while the amorphous systems present inter- 
faces where several types of oxidations and coordination 
levels coexist, the crystalline case presents a situation in 
which all the interface O atoms are single-bonded with 
the Si-NC atoms. Specifically, the Si-atoms at the inter- 





Si32/Si02 


Si 32 -OH 


Si 32 -H 


Crystalline 


2.62 


2.15 


2.78 


Amorphous 


1.19 


0.99 


1.59 



TABLE II: DFT HOMO-LUMO gap values (in eV) for the 
crystalline and amorphous silica, embedded, OH-terminated, 
H-terminated Si32 NC. 
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face of the Siio/Si02 and Si32/Si02 systems are respec- 
tively bonded with 1.6 and 2.8 oxygens in average. It is 
worth to note that a recent size-dependent experimental 
study of Si 2p core-level shift at Si-NC/SiC>2 interface 
showed that the shell region around the Si-NC bordered 
by SiC>2 consists of the three Si suboxide states, Si 1+ , 
Si 2+ , and Si 3+ , whose densities strongly depend on the 
NC size 40 . Generally, such oxydation degree is respon- 
sible for two, competitive, major effects on Eq: while 
from one side a higher oxydation increases Eq , from the 
other side it produces a higher strain on the Si-NC atoms 
that results in a reduction of Eq. Briefly, what emerges 
is that these two contributions are approximately of the 
same order of magnitude, and the dominance of one over 
the other is strictly dependent on the considered system 
(a publication on this topic is currently under develop- 
ment ) . 

In our case, after performing a full relaxation, in or- 
der to remove the strain, of the Siio-OH and Si32-OH 
structures we still obtain an Eq difference of about 1 cV. 
Thus, the increased Eq of the larger system can be re- 
ferred to the greater amount of oxydation of the interface 
silicon atoms. 

The analysis of the valence and conduction band edge 
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FIG. 10: (color online) DFT-RPA imaginary part of the di- 
electric function for the crystalline (dotted) and amorphous 
(solid) Siio (upper) and Si32 (lower) embedded NCs. The 
insets show an enlargement of the spectra at low energies. 



states of the Si32/Si02 system, points out that the for- 
mer are essentially localized in a well defined part of the 
NC, near the interface zone, while the latter are localized 
on the whole volume (from the center to the interface) of 
the embedded structure. 

The effects induced by the amorphization process on 
the optical properties of both Siio/Si02 and Si32/Si02 
systems are depicted in Fig. 10, where the calculated ab- 
sorption spectra for both crystalline (blue line) and amor- 
phous (red line) phases are reported. Here we note that, 
while in the smaller Siio/Si02 system the amorphization 
produces a net change in the shape of the spectrum below 
6 eV, inducing also new strong structures in the energy 
windows between 2 and 4 cV, for the bigger Si32/Si02 
it seems to induce only a simple red-shift. These facts 
enforce our previous hypothesis that for bigger NCs the 
particular configuration of the interface becomes less im- 
portant. This is especially true for the amorphous sys- 
tems, where the diversity of the atomic configurations 
tend to average out the final properties. In order to com- 
pare theoretical and experimental results for the energy 
gap dependency, we introduce a new quantity, E G , which 
is defined as the optical absorption threshold when all the 
transitions with an intensity lower than 1% of the high- 
est peak are neglected. In this way we introduce a sort 
of " instrument resolution" by neglecting very low oscilla- 
tor strength optical transitions (for instance the HOMO- 
LUMO transition in the Siio/Si02 crystalline system). 
The obtained DFT E G values for the embedded systems 
are depicted in Fig. 11. Here E G has been calculated 
for Si-NC containing 2, 5, 10 and 32 atoms. It is clear 
that, while the increasing of the NC size always results in 
a reduction of E G , the effect of the amorphization is to 




dot size (N atoms) 

FIG. 11: (color online) E G versus the Si-NC size and amor- 
phization (embedded systems). N=0 corresponds to the pure 
silica, N= 32, 10, 5, 2 correspond to an average diameters of, 
1.0, 0.7, 0.45, 0.24 nm, respectively. The lines are drawn to 
guide the eye. 
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introduce an additional red-shift that becomes relevant 
for large NCs. For large number of atoms N we expect 
however a vanishing of the red-shift due to the fact that, 
for semiconductors, the main optical properties should 
only depend on the short range-order . 



IV. EXCITONIC AND MANY-BODY 
CORRECTIONS 

Self-Energy and excitonic effects are known to play 
a very important role both in low dimensional systems 
(as the quantum dots) and in 3-dimensional systems 
as the SiC>2. In this last case it is known that DFT 
underestimates the electronic gap of about 5 eV, and 
important excitonic effects arc responsible of the strong 
absorption peak at about 10 eV. The inclusion of the 
many-body effects is thus of fundamental importance 
in order to obtain a better description of the optical 
properties of NCs embedded in SiC>2. Quasi-particle 
effects within the GW approach 42 and excitonic effects 
within the BS equation 3 have been considered in the 
calculation of the Siio/Si02 electronic properties. 
Table III shows the calculated E G within DFT, GW, 
GW+BSE+LF approximations. It is possible to observe 
that, in the ordered (amorphous) case, the inclusion of 
the GW corrections opens up the gap by about 1.9 (1.7) 
eV, while the excitonic and LF correction reduces it by 
about 1.5 (1.6) eV. Thus, the total correction to the 
LDA E G results to be of the order of 0.4 (0.1) eV. We 
note that the final Eg of the crystalline and amorphous 
embedded NCs are quite different, i.e. 2.17 and 1.51 eV, 
respectively. 

The difference between the GW electronic gap and the 
GW+BSE+LF optical excitonic gap gives the exciton 
binding energy Et. Our calculated exciton binding 
energies are quite large: 1.5 eV (crystalline) and 1.6 
eV (amorphous). They are very large if compared with 
that of bulk Si0 2 (almost eV) 44 ~ 46 , bulk Si (~ 15 
meV) or with carbon nanotubes 47,48 where Eb ~ 1 eV, 
but they are similar to those calculated for undoped and 
doped Si-NC 49 ' 50 of similar size and for Si and Ge small 
nanowires 51 ' 52 . 

Fig. 12 shows the calculated DFT-RPA, 

GW+BSE+NLF, and GW+BSE+LF absorption 
spectra for the embedded Siio and a-Siio NCs. The 
results show that the inclusion of the many-body effects 
does not substantially modify the absorption spectra. In 
both cases the energy position of the absorption onset 





DFT 


GW 


GW+BSE+LF 


Crystalline 


1.77 


3.67 


2.17 


Amorphous 


1.41 


3.11 


1.51 



TABLE III: Many-body effects on the gap values (in eV) for 
the crystalline and amorphous embedded Siio dots. 




2 4 6 8 10 



Energy (eV) 

FIG. 12: (color online) DFT-RPA, GW+BSE+NLF (that is, 
without local fields), and GW+BSE+LF (with local fields) 
calculated imaginary part of the dielectric function for crys- 
talline (top) and amorphous (bottom) case. 



is practically not modified (sec insets). Delerue and 
coworkers 53 found that, for Si-NC larger than 1.2 nm, 
the self-energy and Coulomb corrections almost cancel 
each other. In our case (about 0.7nm), this cancellation 
does not completely occur, especially for the crystalline 
system, where self-energy effects dominate. 
Finally, we note that the inclusion of LF contributions 
has dramatic effects on the spectrum of the crystalline 
system, especially in the region between 4 and 6 eV, 
while the amorphous system seems to be affected by a 
simple blue-shift of the absorption threshold. 



V. CONCLUSIONS 

Considering as the /3-cristobalitc phase of the Si02 , like 
as the amorphous phase obtained by classical molecular- 
dynamics, we have generated silicon nanoclusters of dif- 
ferent size and shape, that are embedded in a silica ma- 
trix. The final structures, that have been relaxed with 
the SIESTA method 19 , shows that the surrounding ma- 
trix always produces some strain on the nanocluster, es- 
pecially at the Si/Si02 interface. Then we have compared 



S 



the electronic and optical properties of the embedded 
structures, calculated within the DFT-LDA, with that 
of the free H- and OH- terminated clusters, and that of 
the pure silica. What emerges is that the electronic and 
optical behavior is strongly influenced by the structural 
properties of the considered systems; actually, apart for 
the quantum confinement effect on the energy gap due to 
the different NCs size, the amorphization of the embed- 
ded NC is responsible for the main electronic and optical 
changes of the whole system. 

For the smaller cases, we have also investigated the local 
field, self energy, and binding energy effects. Concern- 
ing the optical properties we have shown that, while the 



many body approach provides only minor corrections, the 
inclusion of local fields is crucial in the analysis of highly 
symmetric systems (crystalline case) but it is of minor 
importance in disordered structures (amorphous case). 
This work shows that it could be worthy to experimen- 
tally investigate the optical properties of the matrix em- 
bedded amorphous nanoclusters 54,55 , that show different 
properties with respect to their crystalline counterparts. 
This work is supported by PRIN2007 and CNR Italia- 
Turchia. We acknowledge CINECA CPU time granted 
by INFM (Progetto Calcolo Parallelo). We thanks A. 
Pedone for the contribution with the modelling of the 
silica glasses. 
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